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Abstract. We investigate the dynamics of a single tracer exploring a course of fixed obstacles 
in the vicinity of the percolation transition for particles confined to the infinite cluster. The 
mean-square displacement displays anomalous transport, which extends to infinite times 
precisely at the critical obstacle density. The slowing down of the diffusion coefficient exhibits 
power-law behavior for densities close to the critical point and we show that the mean-square 
displacement fulfills a scaling hypothesis. Furthermore, we calculate the dynamic conductivity 
as response to an alternating electric field. Last, we discuss the non-gaussian parameter as an 
indicator for heterogeneous dynamics. 
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1. Introduction 

Transport in dense liquids slows down drastically as the packing fraction is increased or the 
temperature is lowered. Approaching the glass transition, transport essentially ceases since 
particles are trapped in long-living cages formed by their neighbors. Then stress fluctuations 
require very long times to relax with the consequence that the viscosity increases drastically. 
The structural arrest of the liquid dynamics is accompanied by a series of characteristic 
features in the correlation functions, such as the emergence of stretched relaxation as terminal 
process and a mesoscopic time window where a critical relaxation appears. A microscopic 
theory for the glass transition has been achieved in terms of the mode-coupling theory [1], 
which has been tested successfully in the last 25 years. 

To further elucidate the microscopic mechanism of the glass transition significant effort 
has been made to confine the liquid to small volumes [2] and thereby to identify the relevant 
length scale of correlated motion. Computer simulations have revealed that walls can 
considerably alter the transport properties [3] and depending on the fluid-wall interactions 
the confinement either facilitates the glass formation [4] or suppresses the tendency towards 
structural arrest [5]. Experiments on colloidal suspensions confined by two quasi-parallel 
walls found an oscillatory dependence of the mean-square displacement perpendicular to the 
walls as the plate distance is varied [6]. A similar variation of the glass transition line due to 
layering has been predicted recently by a mode-coupling theory for confined geometries [7, 8]. 

If the liquid is confined by a random medium, for example a porous host structure, an 
intricate interplay of the freezing of the structural relaxation and a localization transition 
emerges. A series of non-trivial predictions have been made by a generalization of the 
mode-coupling theory accounting for the interaction with frozen obstacles [9-11]. The 
phase diagram for mixtures of mobile and frozen components has been explored by 
computer simulations and semi-quantitative agreement with mode-coupling theory has been 
reported [12-14]. If the density of the mobile component becomes very small, the interaction 
among the mobile particles can be ignored and the dynamics is governed solely by the 
excluded volume due to the frozen matrix. This limiting case is known as the Lorentz 
model which displays a localization transition upon crossing a critical obstacle density. A 
self-consistent mode-coupling kinetic theory [15-17] provides a microscopic description of 
the phase-space correlation functions with no free parameters, which captures quantitatively 
the behavior for low and moderate obstacle densities. In particular, it reproduces the known 
long-time anomalies in the velocity autocorrelation functions [18]. Furthermore, it correctly 
predicts the emergence of anomalous transport at the critical point and the emergence of 
scaling laws. Yet, in the immediate vicinity of the critical point deviations become manifest 
and it appears that mode-coupling theory in its current form is not suited to deal with divergent 
length scales [1]. 

In fact, the localization transition in the Lorentz model is driven by an underlying 
geometric transition of the void space, i.e. the volume accessible to the mobile agents. 
As the obstacle density is increased, the excluded regions grow and start to merge. Then 
configuration space falls apart and a hierarchy of finite pockets coexists with a single spanning 
cluster, which percolates through the entire system [19]. At a certain critical density the 
infinite cluster ceases to have extensive weight rather it becomes a fractal. Then long- 
range transport occurs only on this self-similar structure leading to anomalous transport for 
arbitrarily long times. For even larger densities, all clusters are finite implying that the mean- 
square displacements of individual particles saturate at a finite value. In the vicinity of the 
critical density dynamics scaling is expected to hold [20] with critical exponents from the 
universality class of random resistor networks [21, 22]. 
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Early computer simulations on the ballistic Lorentz model have focused on the non- 
analytic density dependence of the diffusion coefficient [23, 24]. The long-time anomaly in the 
velocity autocorrelation function has been readily observed for moderate scatterer densities 
in two dimensions [23, 24] but the prefactor appeared to be significantly different [25, 26] 
from the theoretical prediction [18]. Only recently, the controversy could be resolved and 
has been shown that the localization transition strongly interferes with the long-time tail [27]. 
For Brownian tracers the velocity autocorrelation function, defined via suitable derivatives 
of the mean-square displacements, has been calculated analytically to first order in the 
scattering density [28] and computer simulations show that the range of validity extends even 
to moderate densities [28, 29]. 

The critical properties of the dynamics of the Lorentz model close to the localization 
transition have been tested only recently, both in three dimensions [30-32] as well as for 
the planar case [29, 33]. Nice agreement with the theoretical predictions has been achieved, 
provided the universal, leading corrections to scaling [34] are taken into account. The scenario 
appears to be relevant also for size-disparate soft-sphere mixtures [35-37] as well as for 
Yukawa mixtures [38] or fast ion diffusion in silica melts [39^-1]. 

The dynamics of the Lorentz model has been characterized so far mainly in terms of all- 
cluster averages, i.e. a tracer is placed at random in the void space and each pocket contributes 
to the mean with a weight given by its volume. However, long-range transport occurs only on 
the infinite cluster below the critical density. Here we elucidate further the dynamics of the 
Lorentz model, Section 2, and study the motion of tracers that are located at the infinite cluster. 
The simulation details are presented in Sec. 3. Then we study the mean-square displacements 
for different obstacle densities (Section 4) and discuss the variation of the corresponding 
diffusivities (Sec. 5). The scaling behavior of the mean-square displacement is elaborated 
in Section 6 where we also discuss the role of the leading corrections. A complementary view 
is to study the frequency-dependent conductivity, which is presented in Sec. 7. The transport 
properties beyond the second moment of the displacement is discussed in Sec. 8 exemplified 
for the non-gaussian parameter. 

2. Lorentz model and continuum percolation 

A minimal model for transport in disordered media has been introduced by Lorentz [42] more 
than a century ago. The idea is that the main contribution to electric resistance arises from the 
interaction of mobile agents (electrons) with frozen impurities in the medium, whereas the 
interaction among the mobile particles is negligible. In its simplest variant these obstacles 
are distributed randomly and independently in the sample and their density n = N/V is 
sufficient to characterize the structure. The mobile agent interacts with the impurities via 
some potential, typically taken as hard-sphere exclusion of radius cr. Then the relevant 
dimensionless control parameter is the reduced obstacle density n* = no 2 . Since the obstacles 
are distributed independently, the regions of excluded volume may overlap and build large 
clusters as the obstacle density increases. Therefore the packing fraction ip, i.e. the volume 
that is inaccessible to tracers, is not merely proportional to the reduced obstacle density n*. 
The relation between (p and n* follows from the following simple argument. The probability 
that a particular obstacle does not occupy a specified point of the large volume V is just 
1 - Ano 3 /2>V. Then the probability that no obstacle covers this point is (1 - Ancr 3 /3V) N , and 
in the limit of large system sizes N — > oo, V — > oo with fixed obstacle density n = N/V, one 
obtains for the packing fraction of the excluded volume 




(1) 
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Only at low obstacle densities <p « Am* 1 3 the probability of mutual overlap of the excluded 
volume can be ignored. For moderate n* the void space consists of a set of isolated bounded 
pockets coexisting with a connected region that spans through the entire system. The latter is 
referred to as the infinite or percolating cluster. 

Above a critical obstacle density n* the infinite cluster ceases to exist and only the 
hierarchy of disconnected finite clusters prevails. The transition is of purely geometric 
origin and is accompanied by a series of power laws reminiscent of a continuous phase 
transition [19]. Directly at the critical density n* c the infinite cluster becomes a self-similar 
structure with fractal dimension^ df = 2.53 implying that its total weight is subextensive. 
The vanishing of its weight as the critical density is approached from below, n* f n*, is 
again determined by a power law ~ {—ef, where e — (n* - n*)/n* denotes the separation 
parameter, a convenient measure for the distance to the critical point. The critical exponent 
(3 assumes the value (3 = 0.41 for three-dimensional percolation. The infinite cluster remains 
fractal below a divergent length scale £ ~ |e|~ v known as the correlation length and appears 
homogeneous only on larger scales. A simple scaling argument reveals that the new critical 
exponent v is connected to the two introduced above by a hyperscaling relation df = d - fi/v 
and in the three-dimensional case one obtains v = 0.88. The structural properties of the 
infinite cluster and the coexisting hierarchy of finite clusters in the immediate vicinity of the 
geometric transition is then characterized by two independent universal critical exponents. 
By universality, introducing short-range correlations in the distribution of the obstacles is 
expected not to alter the nature of the percolation and to leave the values of the exponents 
unchanged. Similarly the scaling functions are considered to be universal and details should 
only enter the amplitudes, e.g. how the correlation length is connected to the radius cr of the 
obstacles. 

The dynamics of particles on such fractal objects becomes self-similar again and 
introduces new critical exponents, which are independent of the ones discussed above. Here, 
we consider a ballistic motion for the tracer scattering elastically at the frozen obstacles. Then 
the magnitude of the velocity v = |v| is conserved and only the direction of v is changed each 
time a scattering event occurs. For a particle initially located on the infinite cluster the linear 
extension R of the probability cloud to find the particle at some distance after an elapsed time 
t is expected to grow as a power law R ~ t 1/dvl directly at the percolation transition. Here the 
walk dimension d w generalizes the concept of Brownian motion where R ~ ? 1/7 also leads 
to a self-similar increase. These arguments suggest that the mean-square displacement of the 
tracer <5r^(f) becomes subdiffusive Sr^t) ~ t 2/d " for long times never reaching the regime 
of diffusion. The mechanism for this phenomenon is encoded in the fractal structure of the 
configuration space and naturally explains the emergence of anomalous transport in complex 
structures. 

The transport properties in the Lorentz model close to the transition point again constitute 
a critical phenomenon with a universality class that is shared by certain random resistor 
networks. There, the fundamental quantity of interest is the vanishing of the conductivity, 
e.g. by cutting randomly bonds of a regular lattice of nodes connected by resistors to their 
nearest neighbors. Then at the percolation transition the residual network becomes a fractal 
with zero linear conductivity. Since the geometric properties of lattice percolation and the 
continuum percolation underlying the localization transition of the Lorentz model have been 
shown to be identical [43], it is suggestive that also the dynamic critical properties are 
related. There is an interesting subtlety though that requires further insight. In continuum 
percolation the obstacles can come arbitrarily close to each other, leading to narrow channels 

$ Recent values for critical exponents for percolation have been collected in Ref. 34 
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through which particles have to squeeze through. These bottlenecks arising due to purely 
geometric reasons weakly connect different compartments of a cluster and the distribution of 
exit probabilities from such compartments is again power-law distributed. Then the mapping 
of the Lorentz model to a random resistor network on a lattice requires to consider power- 
law distributed conductances [21, 22] between nearest neighbors that are diluted until the 
percolation transition is reached. 

3. Simulation details 

Our computer simulations are event-driven molecular dynamics simulations of a single tracer 
particle in a background of randomly distributed, overlapping hard spheres. We use a cubic 
container of linear box size L = lOOcr with periodic boundaries, thus our systems contain 
300,000 to 840,000 obstacles. To study the influence of the finite system size, a number of 
systems with L = 50<x and L = 150cr were also tested. The remaining infinite cluster of void 
space is identified by calculating the Voronoi diagram § of all obstacle centers and analyzing 
the network of edges not overlapping with the spheres. A cluster is regarded as infinite if there 
is an infinite path for a point-like tracer, i.e. the network wraps around of at least one of the 
boundaries [44]. 

For most densities, 16 different obstacle realizations were analyzed with 5 trajectories 
each starting at different points, making a total of 80 independent trajectories per density. 
At n* = 0.83 the number of systems was increased to 80, of which 76 contained an infinite 
cluster, and at n* = 0.84, there was an infinite cluster in 27 out of 100 tested systems. While 
the calculation of trajectories is purely limited by CPU time (one trajectory took about 3 hours 
on a 3.0 GHz Intel Xeon 5160, independently of the system size), the amount of available 
memory limits the number of obstacles for the analysis of the cluster structure (a system of 
1,000,000 obstacles required approximately 3 GB of RAM). 

4. Mean-square displacement 

To characterize the motion of the tracer on the infinite cluster, we have first calculated 
the mean-square displacements 6r^(t) = ([R(f) - R(0)] 2 )oo for obstacle densities up to the 
percolation transition. Here R(f) denotes the position of the tracer and brackets (•)<,<, denote a 
moving time average, an average over different trajectories, as well as the different realizations 
of the disorder. The subscript oo indicates that the particle starts and remains on the infinite 
cluster. 

The results for dr^it) are displayed in Fig. 1 on double logarithmic scales. The mean- 
square displacements extend over more than eight orders of magnitude in time and span 
length scales of a factor of 150 particle diameters thereby displaying essentially no noise. 
For short times the data for all densities collapse on a single straight line which reflects the 
free ballistic motion, 6r^(t) = v 2 t 2 . The mean free path length before the particle encounters 
an obstacle can be estimated by l/ncr 2 implying a collision rate nvcr. Since here we consider 
moderate densities, we measure times in units of to :- cr/v which then is comparable to the 
time of free flight until the next collision event. For rather low scatterer density the mean- 
square displacements cross over directly to a diffusive behavior Sr^if) = bD^t. Increasing 
the obstacle density, the long-time behavior gradually shifts to the right, implying that the 
diffusion coefficient D m is suppressed. 

§ We have employed the free voro++ package initially developed by Chris Rycroft, see http : //math . lbl . gov/ 
voro++/. 
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Time f/v 'cr 

Figure 1. Mean-squared displacements for ballistic tracers on the infinite cluster. Obstacle 
density increases from top to bottom. The straight line corresponds to a power law t 2 ^ w and 
serves as a guide to the eye. 



At a critical density n* = 0.839 the diffusive regime is never reached in the window 
of the simulation time. Rather, the mean-square displacement increases subdiffusively, 
Sr^Jf) ~ f 2 ^ w , over more than six decades in time. The exponent nicely corroborates the 
value dst, = 4.81 which is expected from scaling predictions, compare e.g. [34, 45]. The 
critical density n* also coincides with the point above which no percolating cluster exists for 
large system sizes. 

Just below the critical density the mean-square displacements follow the critical law up 
to a certain time scale and a crossover to diffusive behavior is observed. As the critical density 
is approached the crossover shifts to later times or equivalently to larger distances. The shape 
of the crossover approaches a universal form and scaling behavior is expected, which we 
analyze in Sec. 6. 



5. Diffusion coefficients 



A different representation of the transport properties encoded in the mean-square 
displacements is given in terms of the time-dependent diffusion coefficient 

D ^)=^-Srl{t), (2) 
2d at 

where d — 3 is the dimension of the surrounding space. Figure 2 displays Doa{t) obtained 
by numerical derivatives of the mean-square displacements. Due to the blocking algorithm 
the data are stored essentially on a logarithmic grid such that the derivatives do not introduce 
significant new noise. The short-time ballistic motion becomes manifest as linear increase of 
D M (f) (not shown) before a maximum value is reached. At longer times the time-dependent 
diffusion coefficient drops again, implying that the obstacles induce anticorrelations in the 
velocity autocorrelations. For long times D M (f) saturates and assumes the long-time diffusion 
coefficient D M := D M (f — > oo). The time-dependent diffusion coefficient, defined in Eq. (2), 
converges much more rapidly to its long-time limit than (5r^(?)/6f, since it only probes 
properties at late times. In our simulations, the diffusion coefficient drops by five orders 
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10° 10 2 10 4 10 6 10 8 

Time t/v~ l cr 

Figure 2. Time-dependent diffusion coefficient D m (t) = (l/6)d(5r^,(i)/d; for particles confined 
to the infinite cluster. Obstacle density increases from top to bottom. At the localization 
threshold a power law ( 2 /''»~' becomes manifest, indicated by a straight line as guide to the 
eye. 



of magnitude as the critical obstacle density is approached. Then the subdiffusive motion 
becomes apparent as a power-law decay of Doa(t). 

The extracted diffusion coefficients are expected to vanish as a power law D M (n) ~ (n* - 
n*y>° with a suitable critical exponent p^. Recently, it was shown that the diffusion coefficient 
for an all-cluster average vanishes as D(n) ~ (n* - n*Y with a value p = 2.88 [30, 31]. 
Since only particles that are located in the percolating clusters can contribute to diffusion, 
we anticipate D(n) ~ Dooi^Pooiri) and infer the scaling relation /i^ = p - (3 = 2.47. The 
prediction can be tested in the most sensitive way by rectification plots, see Fig. 3. There 
both D(n) 1 ^ and D co (n) 1 ^ co become straight lines as a function of the scatterer density n*, at 
least for densities above n* = 0.78. Furthermore, they extroplate to the same critical density 
n* = 0.839 which coincides with the density where no cluster percolates through the entire 
system. The scaling behavior is demonstrated also in the inset of Fig. 3 where we show that 
D(n)\e\^' and D co (n)|e| _ ' J " approach constants as the separation parameter e becomes smaller. 



6. Scaling behavior 

The similarity of the mean-square displacements and the power-law dependence of the 
diffusion coefficient in the vicinity of the percolation transition suggests a scaling law. Then 
6r^(t; n) is no longer a function of both time and density but the density dependence can be 
hidden by suitable shifts in time and length scales. Directly at n* the critical law dr^ ~ t 2ld ~" 
persists for all times suggesting the scaling hypothesis 

6rl(t\e) = ? ld "6rl{tlt x ), (3) 

where the time scale t x diverges for n* f n*. For short rescaled times, the scaling function 
assumes a constant Sf^it <§c 1) = const., since the particle has not explored the infinite cluster 
on large enough scales to sense deviations from a fractal structure. On the other hand, for large 
rescaled time the mean-square displacement increases linearly implying <Sr^(f » 1) ~ f 1-2 '"". 
Furthermore, one expects that the crossover from subdiffusive to linear increase of the 
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Figure 3. Rectification of the diffusion coefficients for the infinite cluster and the all- 
cluster average D. Both diffusion coefficients vanish at the same critical obstacle density 
n* = 0.839. Inset: Scaling behavior of the diffusion coefficients with the separation parameter. 



mean-square displacement dr^if) occurs on a length scale given by the correlation length 
drlaitx) ~ £ 2 , suggesting a power-law divergent crossover time scale t x ~ if" ~ (-e)~ vrfw . For 
the long-time diffusion coefficient one infers D M ~ t 2 J d ^^ 1 ~ £ 2_d » and we obtain a relation 
between the critical exponents p^ = v(c/ w - 2) [45]. The values that we have extracted from 
the subdiffusive increase of the mean-square displacement (5r^(f) and the vanishing of the 
corresponding diffusion coefficient D M nicely corroborate the exponent relation. 

Earlier it has been demonstrated that the all-cluster averaged mean-square displacement 
also displays subdiffusive behavior 6r 2 (f) ~ t 2 ^~~ at the critical obstacle density [30, 31], yet 
with a different dynamic exponent z = 6.25. Here the self-similar hierarchy of finite clusters 
contributes to the mean-square displacement. Below the critical density n* < n* a regime 
is reached where 5r 2 (t) increases again linearly, 6r 2 {t) - 6Dt, and D ~ (-eY,p — fi x + ft 
decreases faster than D M ~ (—eY°° since only the particles initially on the infinite cluster 
contribute to long-range transport. The time needed for these particles to explore the 
infinite cluster is still given by t x ~ f w and and the scaling law 5r 2 {t;e) = t 2 ^ z 5r 2 (t/t x ) 
is anticipated. The scaling function 6r(i) encodes the crossover from anomalous transport 
5f 2 (i <sc 1) = const, to the linear regime 6f 2 {i » 1) ~ f 1-2 ^. In particular, the scaling 
hypothesis predicts that the diffusion coefficient vanishes as D ~ t x ~ {—ef which leads 
to the exponent relation p. = vc/ w (l - 2/z). Combining with the scaling relation obtained so 
far, one finds for the dynamic exponent z = vd w /(y - ft/2) and the measured value z = 6.25 
is in excellent agreement with the prediction. Second, the length scale t where the crossover 
from subdiffusive to linear increase occurs, I 2 ~ 6r 2 (t x ) is given by I ~ t l J z . Employing the 
scaling relations derived so far, one finds I ~ (-eY v+l312 and it is determined entirely by critical 
exponents for structural properties of the percolation transition. In particular it diverges slower 
than the correlation length £ ~ (-e)~ v and in fact constitutes the second moment of the size 
distribution of the clusters [19] and the measured values of £ were shown to corroborate the 
scaling hypothesis [30, 31]. 

The consistency of the values of the exponents expected from the exponent relations is 
encouraging to test the scaling hypothesis, Eq. (3) itself. Thus upon dividing by the critical 
law ? 2 / rfw and measuring times in units of the crossover time t x ~ f" all data sets are expected 




Figure 4. Scaling behavior of the mean-squared displacement on the infinite cluster. Left 
panel: Upon rescaling with the crossover time t x '■= to(-£)~ vclw with to;= cr/v the data 
converge to a master curve as the critical density is approached. Right panel: Leading 
corrections are approximately taken into account by adjusting a single nonuniversal amplitude 



to collapse onto a single master curve Sr^it), see Fig. 4. The curves all display the same long- 
time asymptote and exhibit a satisfactory data collapse. However, this reflects only that for all 
data sets our simulations reach the diffusive regime and that the diffusion coefficients follow 
the scaling prediction. Deviations from scaling become manifest at low rescaled times and 
here one infers only that the data sets approach a master curve as the separation parameter 
to the critical point is decreased. The crossing of the asymptote occurs where the scaling 
regime merges continuously with the microscopic time scales. The fanning out of the data 
for short rescaled times has also been observed for the all-cluster averaged mean-square 
displacements [30, 31] and it was argued that corrections to scaling should be taken into 
account. 

The idea in the renormalization group is that rescaling of length scales corresponds to 
changing the coupling constants in some large dimensional space of parameters. The critical 
point corresponds to a fixed point of the renormalization group flow and the properties of 
this fixed point determine the values of the critical exponents. The scaling functions encode 
how the critical point repels trajectories that are close to but not on the critical manifold. 
In particular, the corrections to scaling are universal again and contain the information how 
the scaling regime is approached or how microscopic details become irrelevant. These ideas 
suggest that the scaling hypothesis, Eq. 3 can be extended to account for generic correction 
terms [34] 

Sr\t- e) = t 2/d -6ri(t) [1 + r v A 0O (f)] , f = t/t x . (4) 

The exponent y is again universal and has been determined to y — 0.34 for the three- 
dimensional Lorentz model from a similar scaling hypothesis for the all-cluster averaged 
mean-square displacement [30, 31]. The new scaling function assumes a constant for small 
rescaled times C M := A M (f <K 1) and behaves like a power-law at large scaling times 
AooCf » 1) = P. In principle, the amplitude can be inferred from the mean-square 
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displacement at the critical obstacle density 

<Jr 2 (f;e = 0) ~ f 2A/ " [1 +Coor- 



(5) 



Similarly, the amplitude encodes the corrections to scaling for the diffusion coefficient 



Since the detailed shape of Aoo(f) is unknown, one has to rely on simple approximation 
schemes. A convenient interpolation function could be A 00 (f) = C m + C^P . Since in our 
data the corrections to scaling are most significant at short rescaled times, we ignore the 
corrections to the diffusion coefficient C^°' « 0, and approximate A M (f) = Coo for all times 
i. The remaining unknown coefficient Coo is adjusted to yield an optimal data collapse, see 
Fig. 4. Choosing Coo = -0.53 ? Q the data nicely collapse on a single master curve although the 
improvement is less satisfactory than for the all-cluster average [30, 31] or in two-dimensional 
transport on a lattice [34]. The reason may be traced back to the slower convergence of the 
diffusion coefficients identified in the inset of Fig. 3. 

7. Conductivity 

In applications, such as ion transport in disordered solids, the basic quantity of interest 
is the frequency-dependent complex conductivity <x(to). Applying a small homogeneous 
alternating electric field E(f) = Re[E(to)e~ lw '] to the sample the resulting electric current 
density will alternate with the same frequency j(f) = Re[j(w)e _lw '] in the linear response and 
the appropriate frequency-dependent generalization of Ohm's law reads j(to) = <x(to)E(to). 
The real part cr'(to) = Re[cr(<y)] > describes the dissipation due to friction, whereas the 
imaginary part cr"(o)) = Im[cr(to)] is connected to the energy stored in the system. 

The current density constitutes a collective property of the sample but for dilute charge 
carrier concentration n lon and charge density qn\ on the interaction between the mobile charges 
may be ignored. Then each ion experiences a frequency-dependent force F(o>) = qE(co) 
resulting in an oscillatory velocity \(co) = /i(to)F(to), where p(a>) denotes the frequency 
dependent mobility, sometimes also referred to as admittance [46^-8]. Since for independent 
ions the total current density is given by j(to) = ni on q\(cj) the frequency-dependent 
conductivity can be related to the frequency-dependent mobility /i(to) by o-(a>) = q 2 ni on p(a>). 
By the well-known fluctuation-dissipation theorem [49] yu(to) = Z(aj)/k B T where 



is the one-sided Fourier transform of the velocity autocorrelation function. In particular, the 
Green-Kubo relation is obtained as low-frequency limit D = Z(cj — > 0) and one recovers the 
Stokes-Einstein-Sutherland relation /i = D/k B T. In summary, Z(co) encodes the frequency- 
dependent response to an alternating field and coincides up to well-understood factors with 
the frequency-dependent conductivity. 

Here we single out the conductivity for the particles located on the infinite cluster only. 
The experimental condition in mind is that the system initially does not contain ions but they 
are inserted in the infinite cluster by wiring the surfaces. Then the measured conductivity is 
essentially Zoo(to), viz. the one-side Fourier transform of the velocity autocorrelation, where 
the average is restricted to particles on the spanning cluster. 

The numerical evaluation of Zoo (to) is achieved upon integrating by parts 



Doo(e) ~ c-er [i + c£V(- e r yvdw ] 



(6) 




(7) 




(8) 
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Figure 5. Real part of the frequency-dependent conductivity for particles restricted to the 
infinite cluster. Obstacle density increases from top to bottom. At the critical obstacle density 
n* the conductivity follows a power law indicated by the black line as guide to the eye. 

and the Fourier transform can be easily performed relying on the simplified Filon 
algorithm [50]. The loss part of the conductivity Re[Zco(w)] is displayed in Fig. 5 for 
a frequency window covering more than five decades and three orders of magnitude in 
amplitude. At the critical density «*, anomalous transport becomes manifest in terms of a 
power-law w I ~ 2/,rf *. This fractal behavior is inherited directly from the power-law decay of 
the time-dependent diffusion coefficient D^if) ~ t 2 ^* 1 ^ 1 or equivalently by the subdiffusive 
increase of dr^it) ~ f 2 / d ». 

For lower scatterer densities n* < n* the curves display a finite low-frequency plateau 
which decreases drastically as the percolation threshold is approached. One identifies this 
low-frequency limit with the diffusion coefficient as = Zo<,(a>) or with the dc conductivity, 
respectively. In the vicinity of n* the curves approach the power law increase w 1-2 ^" in 
an intermediate time window that gradually opens as n* f n*. The crossover from the dc 
conductivity plateau to the critical behavior occurs at a crossover frequency co x that shifts to 
lower frequencies as the density comes closer to the percolation transition. Upon inspection 
of the curves one anticipates scaling behavior 

Z M (w; e) = w 1 - 2/d »Z M (w), Cj = lo/lo x , (9) 

where the scaling function Z.cx,(fi>) encodes the crossover. By construction » 1) = 

const., whereas the behavior for low rescaled frequencies follows from the requirement that 
the conductivity remains finite, Z^iCS) ~ &> 2 ' <iw_1 . Consequently the diffusion coefficient 
scales as D M ~ wjr 2 ^ w and from the previously derived scaling relations one infers a> x ~ 
(-e) vd '\ In particular, the crossover frequency of the conductivity scales as the inverse of the 
crossover frequency in the mean-square displacement co x ~ r x x as anticipated from the Fourier 
transforms. 

8. Non-gaussian parameter 

For ordinary diffusion, the statistics of the increments AR(f) := R(f) - R(0) is given by a 
gaussian distribution, characterized entirely by the first two cumulants. The first vanishes by 
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rotational symmetry (AR(f)) = 0, whereas the second is simply the mean-square displacement 
6r 2 (t) = ([AR(f)] 2 ) = 6Dt. The basic assumption underlying the picture of ordinary diffusion 
is that the increments are independent at least for sufficiently long lag time t and the central 
limit theorem applies. Then by independence 5r 2 (f) ~ t increases linearly and the diffusion 
coefficient D characterizes the variance of the increments for a unit time. 

In the case of subdiffusion the central limit theorem is necessarily violated, and the higher 
cumulants encode valuable information on the statistics of the increments. By rotational 
symmetry all odd moments still vanish identically, and here we focus on the next higher 
nontrivial moment 6r^,(t) := (|AR(f)| 4 ) M , where we restrict again the average to particles 
confined to the infinite cluster. The simulation results are displayed in Fig. 6 and one identifies 
a power-law increase at the critical density <5r^(f) ~ f 4 ^ w over more than six decades in time. 
The data for very long times are presumably affected by statistical fluctuations, viz. whether 
an infinite cluster is present or not. Since we have averaged the data over those realizations that 
in fact display such an infinite cluster, the curves are expected to display some statistical bias. 
For densities below the percolation threshold the data reach long-time asymptotes 6r^(t) ~ t 2 
as expected for diffusive motion, and the crossover time is identified with the one of the mean- 
square displacements (not shown). 

A more sensitive quantity characterizing the differences to a gaussian statistics of the 
increments is the dimensionless non-gaussian parameter 

a t\t) = \^r 2 -h (10) 

where the prefactors are chosen such that for a three-dimensional gaussian transport ct^ (?) = 
for all times. The non-gaussian parameter has been identified in the context of supercooled 
liquids and dense colloidal suspensions as an indicator for heterogeneous dynamics, i.e. a 
growth of the non-gaussian parameter in an intermediate time window suggests that several 
processes with a broad spectrum of individual relaxation times are present. In Fig. 7 we 
display the non-gaussian parameter for particles confined to the infinite cluster as well as the 
one corresponding to the all-cluster-averaged motion az(f) = (3/5)5r 4 (t)/[6r 2 (t)] 2 - 1. The 
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Figure 7. Non-gaussian parameter for particles confined to the percolating cluster. The full 
lines correspond to particles restricted to the percolating cluster, whereas the broken lines are 
the all-cluster-averaged dynamics. 

most prominent difference is that aq (i) approaches zero for long times for densities below 
the localization threshold, whereas ai{t — > oo) saturates at finite values for n* < n*, which 
grow upon approaching the transition. Directly at the critical point a2(t) is predicted to diverge 
as a power law [30], in contrast to the finite limiting value a 2 (t — > oo). The data for aq (f) 
at n* are sensitive to simulation details and appear to increase drastically at late times. The 
inset in Fig. 7 tests the long-time behavior of aq (f) for different box sizes and we conclude 
that the observed increase is indeed a finite-size effect. 

The decay to zero at late times of ai (?) for densities below n*. suggests that the statistics 
of increments, once the dynamics is restricted to the percolating cluster, becomes a gaussian. 
To characterize the motion on the infinite cluster beyond the lowest moments one should 
calculate the van-Hove correlation function Gca(r, t) :- (6[r - AR(f)])co, which represents the 
probability for the tracer to have moved a distance r in a lag time t on the infinite cluster. In 
particular, all moments of the distribution of AR(f) can be calculated. Scaling suggests that at 
criticality 

Goc(r, t) = r~ d G m (?) , t = t/t x ~ (11) 

implying that the even moments scale as 5r 2n (t) ~ f 2 "^" [45] with nontrivial amplitudes 
determined from the scaling function Gco(0- Below n* the infinite cluster appears 
homogeneous on length scales larger than the correlation length ^ and one expects to recover 
the gaussian propagator for the motion. 

Let us emphasize again that the behavior of the all-cluster-averaged motion is 
qualitatively different from the infinite-cluster-only dynamics. The latter displays anomalous 
transport with propagators that deviate from simple gaussians for length scales r <K ^ and 
time scales t <K t x , whereas the long-time large-distance motion is described again by simple 
diffusion as anticipated by the central limit theorem. In contrast, the all-cluster-averaged 
motion never reaches the diffusion law, although the corresponding mean-square displacement 
increases linearly for times larger than the crossover time t x . The frozen disorder implies that 
configuration space is compartmentalized for all densities such that the particles confined 
to some finite cluster remain there for all times. Thus the all-cluster-averaged dynamics is 
intrinsically heterogeneous at all densities and all times. 



REFERENCES 



14 



9. Summary and Conclusion 

Long-range transport in compartementalized solids is restricted to the infinite cluster 
percolating through the entire system. We have established within the Lorentz model as 
minimal model for such heterogeneous media that the arrest of the dynamics coincides with 
the percolation threshold. Then the infinite cluster becomes a fractal object and transport on 
it becomes anomalous. Here we have employed large scale computer simulations to test the 
dynamics at the localization threshold and in its immediate vicinity. We have corroborated 
that the dynamics can be rationalized by a critical phenomenon, similar to a continuous phase 
transition, and tested various aspects of the critical behavior. In particular, we find that three 
independent critical exponents are necessary to obtain a complete description of the dynamics. 
Here the exponent /3 quantifies the vanishing of the weight of the infinite cluster, v the increase 
of the correlation length, and d w characterizes the mean-square displacement at the transition 
point. The emergence of subdiffusive transport in the mean-square displacement and scaling 
behavior is accompanied by a vanishing diffusion coefficient as well as the appearance of 
an anomalous conductivity. Furthermore we have discussed the violation of the central limit 
theorem directly at the critical point and its consequences for the non-gaussian parameter. 

The present study targets ion transport in a quasi-frozen, three-dimensional environment. 
In these systems the disorder certainly contains correlations as the environment consists of 
some frozen-in configuration of a densely-packed strongly-interacting liquid. Thus it would 
be interesting to study transport also in disordered systems displaying some short-range 
correlations and see, whether the general scenario presented here qualitatively applies, in 
particular if the values of the dynamic exponents remain unaffected by the correlations. Since 
it has been shown, that the value of the conductivity exponent is sensitive to the distribution of 
the weak conductances between different compartments of the clusters, the answer depends 
on how the local geometry affects the statistics of these bottlenecks and how they are tested by 
the particle's dynamics. Second, to make a comparison to realistic ion-conducting melts, one 
should take the finite lifetime of the local structures into account. Since the obstacles slowly 
rearrange, once in a while blocked pockets will be connected to the infinite cluster and new 
dynamical processes enter the problem. 

It would be interesting to study also the dynamics on the infinite cluster for two- 
dimensional systems as minimal model for biological plasma membranes. Experiments on 
protein transport in the model cell membranes by fluorescence correlation spectroscopy (FCS) 
have revealed anomalous transport [51-53], which has been interpreted as a manifestation of 
crowding. Computer simulations modeling the obstructed motion in such membranes indeed 
observe anomalous transport [33, 54-58] and have identified the localization transition as a 
generic candidate to explain the experimentally observed findings. 
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